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ABSTRACT 

Selected examples are described in which NASTRAN is used as an analysis 
research tool for composite mechanics and for composite structural components. 
The examples were selected to illustrate the importance of using NASTRAN as an 
analysis tool in this rapidly advancing field. The results obtained demon^ 
strata rather convincingly the versatility and effectiveness of NASTRAN in such 
applications. 


INTRODUCTION 

This paper describes selected examples in which NASTRAN was used as an an- 
alysis tool for composite mechanics and for composite structural components. 
Typical examples for composite mechanics include: stress analysis of the load 

transfer region of off-axis (with respect to fiber direction) composite spec- 
imens, the vibration response of specimens with defects, the effects of in-plane 
and out-of-plane (plane of specimen) eccentricities on test results, the stress 
variation in the load transfer region of thin composite tubular specimens, the 
stress distribution in the near and far field regions of a composite cantilev- 
ered plate subjected to a concentrated load, and the determination of the bend- 
ing modulus in composite specimens with non-uniform thickness. The examples 
for composite structures include; stress analysis of composite blades for air- 
craft turbine engines, and stress and vibration analysis of a composite thin 
shell airfoil section. 

The description of each example includes; the objective of the study, a 
brief background of its genesis, a schematic of the geometry, the finite ele- 
ment representation (nxnaber of grid points, number and type of elements, com- 
puter running times), typical results in graphical and/or tabular form, and 
comparisons with measured data when available. 

OFF-AXIS COMPOSITE TENSILE SPECIMENS 

The objective of using NASTRAN in this study was to predict the stress 
state near the end tabs and at mid length of off-axis composite specimens test- 
ed in tension. The results were subsequently used for comparisons with mea- 
sured data and for identifying possible out-of-plane eccentricities that might 
be present during testing. 

The geometry and Instrumentation of the specimen is shown in figure 1. 

This specimen was made from high modulus graphite-fiber epoxy resin composite 
(Modnor l/epoxy (MOD I/E)). The laminate configuration of the specimen con- 
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sisted of 8-plies COg], all oriented in the same direction. The specimen 
thickness was about .057 inch and it was loaded at 30° to the fiber direction. 

The NASTRAN model is shown in figure 2. The NASTRAN model consisted of 
657 nodes and 576 quadralateral plate (CQUAD2) elements which included the tap- 
ered portion of the reinforcing end tabs. Note that the finite element repre- 
sentation includes two groups of elements. At each end, the elements are 
0.0625 in. long; these represent the tapered portion of the reinforcing tabs and 
the first quarter inch segment of the test section which is the site of the top 
strain gages. The remaining elements of the representation are 0.125 inch in 
length. All elements for this model are 0.0625 wide. The element size was 
made small enough to study the zones where the strain gages were located on the 
actual specimen. The material properties required for NASTRAN were generated 
using the composite mechanics code (ref. 1). Typical CPU (central processing 
unit) time was 18 minutes for the UNIVAC 1110. 

Axial strains across the width of the 30° off-axis specimen as determined 
by NASTRAN are shown in figure 3 along with the measured fracture strains. Mid- 
length measured strains at fracture ranged from 0.365 to 0.380 percent while 
those at the tab end were 0.395 and 0.424 percent. The NASTRAN analysis axial 
strains across the specimen midlength were between 0.374 and 0.426 percent and 
those at the tab-end were between 0.353 and 0.447 percent. A glance of fig- 
ure 3 shows that NASTRAN predicted strains are a little higher than the mea- 
sured strains at the specimen midlength and lower than the measured strains 
near the specimen end. Although not shown in figure 3, NASTRAN strain predic- 
tions for a section of the specimen lying midway between the tab-end and spec- 
imen midlength were between the two curves, as would be expected. The discrep- 
ancy between the strains measured at midlength and tab-end of the specimens are 
not accounted for by the NASTRAN analysis. 

This leads to the conclusion that the discrepancy is probably produced by 
possible out-of-plane load eccentricities which were investigated and are de- 
scribed in the next section. The important point to be noted is the use of 
NASTRAN in identifying possible test difficulties. 

OUT-OF-PLANE ECCENTRICITIES IN COMPOSITE OFF-AXIS SPECIMENS 

The objective of using NASTRAN in this study was to assess whether pos- 
sible out-of-plane of the specimen eccentricities such as bending or twisting 
could produce part or all of the discrepancy described in the previous section. 

The geometry of the specimen and the NASTRAN model are the same as de- 
scribed in the previous section (figs. 1 and 2, respectively). For this study, 
the specimen was first loaded with an out-of-plane (plane of specimen) bending 
moment and then with a twisting moment at one end of the specimen. The ends of 
the specimen were constrained to lie in the same plane for the out-of-plane 
bending case. The center line of the specimen was constrained to lie in the 
same plane for the twisting case. 

The out-of-plane bending effects predicted by NASTRAN on axial strain for 
the 30 off-axis specimen subjected to a 100 inch-pound bending moment are 
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shown in figure 4. Those for torsion are shown in figure 5 (100 in. -lb twist- 
ing moment) . 

At specimen midlength, the twisting effects are very small, figure 5, but 
near the grip they are important and are of opposite signs on the two edges of 
the specimen. The out-of-plane bending effects shown in figure 4 are moderate 
but all in the same direction at specimen midlength. Near the tab they are of 
the opposite sign and are much larger at one edge of the specimen than at the 
other. The combined effect of twisting and out-of-plane bending near the tab 
end of the tensile specimen could be sizeable and data provided by a strain 
gage at such a location could be erroneously interpreted. 

The important conclusion from the previous discussion is; out-of-plane 
eccentricities of tensile specimens could give rise to very large strains near 
the end tabs of the specimen. This may be a cause of the frequent failure of 
this type specimen near the tab. It is significant to note that to have ar- 
rived at the same interpretation through the use of strain gages would have re- 
quired a large number of strain gages which would have been impractical. 
Therefore, the effectiveness of using NASTRAN in this study as a research tool 
is obvious. 


VIBRATION OF COMPOSITE SPECIMENS WITH DEFECTS 

The objective of using NASTRAN in this study was to predict the effects 
of defects in composite specimens on the free vibration modes of these speci- 
mens. For this purpose, the vibration modes of composite specimens with pro- 
gressively larger through-the-thickness defects were determined. 

The specimen geometry is shown in figure 6. The specimen analyzed was 
assumed to be made from Thornell 75/Epoxy (T75/E)[±45, 90, 90, 0, 0]g compos- 
ite. The defects considered were slits 1/4, 1/2, and 3/4 of the specimen width, 
in length, respectively, and 1/4 inch wide. 

The NASTRAN model is shown in figure 7. The model consists of 297 nodes 
and 256 CQUAD2 elements. The through-the-thickness defect was simulated by re- 
moving elements 124 and 125 for the 1/4, 123 to 126 for the 1/2, and 122 to 127 
for the 3/4 size defect, respectively. The material properties required for 
NASTRAN were obtained from the composite mechanics computer code (ref. 1). 
Typical CPU times were 17 minutes in the UNIVAC 1110. 

Free vibrational frequencies were calculated for an undamaged specimen and 
for a specimen with defects at midlength which spanned 1/4, 1/2, and 3/4 of the 
width of the specimen by removing the connection cards for the appropriate el- 
ements. Table I shows the frequencies for the first four modes for each con- 
dition. The frequencies are lowered by the defect. Table II shows the percent 
reduction in specimen area resulting from the defects of the three sizes as well 
as the relative changes in the corresponding frequencies. A defect spanning 
one-fourth of the specimen width which removes 0.78 percent of the bulk of the 
material in the specimen reduces the vibrational frequency by no more than 1.2 
percent in the case of the first four m.odes. This would represent a very large 
crack in an actual component, which leads to the conclusion that it would take 
a large defect or Imperfection to affect the integrated stiffness and mass of 


383 


the component. Also the vibration frequencies are not very sensitive to small 
defects in composite structures- Figure 8 shows the fourth mode shapes for an 
undamaged composite specimen and the 1/2 size defect (1.56 percent of the spec- 
imen area). The mode shapes are almost identical and the frequency drop due to 
this large defect is less than 3.0 percent. 

The important conclusion from this study is that the vibration frequency 
measurement used as an NDE procedure for evaluating a composite structure with 
defects of the type investigated is not likely to prove a practical procedure. 
This example illustrates the use of NASTRAN as a research analysis tool for as- 
sessing the practicality of such NDE procedures. 

TRANSITION REGION STRESSES IN TUBULAR COMPOSITE SPECIMENS 

The objective of using NASTRAN in this study was to predict the stress 
state near the grip area (transition region) in thin composite tubular speci- 
mens. The predicted stress state was used to (1) compare it with measured data 
and (2) more importantly, evolve grip designs which minimize the stress in the 
transition region. Thin tubular specimens are important in studying failure 
and failure mechanisms of composites subjected to controlled combinations of 
multlaxial loading, references 2 and 3. 

Numerous grip failures were occurring in tests of tubular unidirectional 
composite test specimens. The specimens were potted into metal grips with an 
epoxy resin as shown in figure 9. One of these specimens was instrumented with 
strain gages in the grip transition region so that predicted results using 
NASTRAN could be compared with measured data. The geometry of the specimen 
which was used in the NASTRAN stress analysis is shown in figure 10. The com- 
posite material used in the analysis was MOD I/E with a laminate configuration 
of [Og]. The material properties required for use in NASTRAN were determined 

from data provided by the material supplier. 

A schematic of the NASTRAN model is shown in figure 11. The model con- 
sisted of 576 nodes and 552 CQUAD2 elements, 24 around the circumference and 
23 along the length. These elements were used to model the tube wall and the 
potting material. The elements representing the potting material were normal to 
those for the tube wall. Typical computer CPU times were 21 minutes in the 
UNIVAC 1106. 

Figure 12 compares the experimental and predicted results. The good 
agreement between the predicted results and the experimental data gave confi- 
dence in the NASTRAN model. It was then used to determine the effect on grip 
transition stresses of different grip designs and potting material properties. 
With the aid of this model it is possible to Identify grip designs with re- 
duced transition stresses and, thereby, evolve a more practical grip design* 

BENDING MODULUS OF COMPOSITE ANGLEPLIED LAMINATES WITH SOME 

THICKNESS VARIATION 

The objective of using NASTRAN in this study was to identify an effective 
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means for determining the bending modulus of specimens from composite angle- 
plied laminates with variations in thickness. The bending modulus and the usual 
elastic modulus, measured in a tensile test, are different for angleplied lam- 
inates in general. 

The material was MOD I/E and was in the form of tensile specimens with end 
tabs already in place and the thickness of the test section varied by as much 
as 7 percent. For cantilever bending, stiffness is related to fundamental fre- 
quency by the following expression: 

E=3.2f2^ (1) 

where E is the bending modulus, f^ is the fundamental frequency, m is the 
mass per unit length, L is the length of the cantilever, and I is the moment 
of inertia. 

The fundamental frequency was determined experimentally. To relate this 
frequency to bending modulus, the specimen was modeled for NASTRAN taking into 
account the thickness variation and the end tab. Assuming that the bending mod- 
ulus was constant along the length of the specimen, the modulus was determined 
from the relationship 



where Eg is the effective bending modulus, is a trial bending modulus used 
In the NASTRAN model, f^ is the experimentally determined fundamental bending 
frequency, and fjj is the fundamental frequency obtained from the NASTRAN mod- 
el. Schematics of the tensile specimen and the NASTRAN model are shown in fig- 
ure 13. Bending moduli determined from the above procedure for four quasi- 
isotropic angleplied laminates with different laminate configurations are shown 
in table III. This example Illustrates the versatility of NASTRAN. 

HIGH VELOCITY IMPACT COMPOSITE CANTILEVER SPECIMENS 

The objective of the NASTRAN stress analysis of this specimen was to de- 
termine high stress regions and stress types under point load. The information 
obtained was used to: (1) guide the selection of sites for placing strain 

gages to measure the high velocity impact response of composite centilevers and 
(2) to help interpret some of the high velocity impact data generated in ref- 
erence 4 under contract to NASA LeRC. 

A photograph of an impacted specimen reported in reference 4 is shown in 
figure 14. The geometry of the specimen is depicted in figure 15. The NASTRAN 
model is shown in figure 16. The NASTRAN model consisted of 333 nodes and 
288 CQUAD2 elements. Typical CPU times were about 4.5 minutes in the UNIVAC 
1110. The specimen, for which results will be presented, was made from intraply 
hybrid composite with the following constituents: 75-percent type AS graphite 
fiber-epoxy resin composite (AS/E) and 25-percent S-glass - epoxy resin com- 


posite (S-Gl/E) . The laminate configuration for the specimen was [±40, 0, 10, 
" 10]2 which is representative of composite blade airfoil laminate config- 
urations to minimize residual stresses. The material properties required for 
use in NASTRAN were generated using the composite mechanics code, ref. 1. 

Predicted stresses on the back surface of the composite specimen due to a 
unit point load at the impact point are shown in stress contour plots as fol- 
lows (refer to fig. 16 for directions): spanwise stress (a ) in figure 17(a), 

chordwise stress (Cy) in figure 17(b), and in-plane shear stress (o^y) in fig-- 
ure 17(c). The important points to be observed from these stress contour plots 
are: 

1. High spanwise stresses at the support near the specimen center, fig- 

ure 17(a) 

2. High chordwise stresses near the specimen center, in-board from the im- 

pact point, figure 17(b) 

3. High in-plane shear stresses at the specimen edges near the support, 

figure 17(c) 

The important conclusion from the above observations is that high stress 
regions for these types of specimens are readily identified using NASTRAN. 
Strain gages should be placed at these regions for obtaining optimum informa- 
tion. It is also important to note that the estimated maximum impact force 
built up during these tests is about 1000 pounds for 700 feet per second pro- 
jectile velocity within 50 microseconds from initial contact (ref. 5). This 
magnitude of impact force will produce high tensile chordwise stresses in the 
back surface of the specimen near the impact point, figure 17(b), which, in 
turn, will produce splitting and delaminations at this region. This is consis- 
tent with the test result exhibited in figure 14. 

HIGH-TIP-SPEED COMPOSITE FAN BLADE 

The original objective of this investigation was to develop a computerized 
capability which couples NASTRAN with composite mechanics for the structural 
and stress analysis of composite fan blades for aircraft engines. A descrip- 
tion of this capability and results obtained therefrom have been reported pre- 
viously in references 6, 7, and 8. Herein, predicted stresses at design loads 
are compared with fracture stresses in tensile specimens cut from the critical- 
ly stressed region in the composite blade. The design loads include aerody- 
namic pressure and temperature, and centrifugal forces resulting from 2200 feet 
per second tip speed with a mean tip radius of 16.4 inches (ref. 6). 

A photograph of the fabricated composite blade is shown in figure 18. The 
blade is twisted, cambered, and tapered (both spanwise and chordwise). The 
blade was made from high-tensile-strength graphite fiber in Kerimid polyimide 
matrix (HTS/K601). The blade consisted (blade laminate configuration) of 73 
plies at its thickest portion 30-percent of which were oriented at ±40^ to the 
radial direction at the surface (shell plies) and 70-percent were oriented at 
0 to the radial direction at the center (core plies). This type of laminate 
configuration is usually called shell/core. In addition, the blade had ±20° 
transition plies between the ±40° shell plies and the 0^ core plies, and also 
±70 plies near the tip for increased flutter resistance. For additional de— 
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tailed description see references 6 and 8. 

The NASTRAN model of the composite blade is shown in figure 19. The model 
consisted of 299 nodes and 531 triangular plate (CTRIA2) elements. Typical CPU 
times were 11.5 minutes in the UNIVAC 1106. The average stresses predicted by 
NASTRAN for the design load are shown in a contour plot in figure 20. Stresses 
from the high stress region points A, B, and C in figure 20 are tabulated in 
table IV. Also shown in this table are the material fracture stresses measured 
using tensile specimens cut from the same region from fabricated blades (in situ 
material fracture stresses). As can be seen in table IV, the predicted stresses 
on the pressure surface are higher than the fracture stress. As a point of in- 
terest, some of these blades failed below or at design load during spin tests. 

The important conclusion from the previous discussion is that NASTRAN 
coupled with composite mechanics can be used to predict fracture stress in struc- 
tural components with complex geometry in anisotropic heterogeneous materials. 

The significant point to be emphasized here is that predicted stresses must be 
compared with in situ material fracture stresses for reliable assessment of the 
design. 


FIBER COMPOSITE THIN SHELL AIRFOIL 

The objective of this study was to select on a preliminary design basis a 
laminate configuration for NACA 64A010 type airfoils for wind tunnel wing sim- 
ulation studies. The design load conditions were: (1) dynamic pressure typical 

for this type of airfoil, (2) first mode natural frequency of 60 Hz or greater, 
and (3) a maximum weight of 20 pounds. The margin of safety on stresses due to 
dynamic pressure was about 5, to account for possible fatigue, and on the first 
mode frequency was about 1.5. 

The geometry of the airfoil section is shown in figure 21. Note the air- 
foil is a thin shell so that the weight is kept to a minimum. The composite 
material selected was graphite/epoxy (AS/E) 32 plies (0.005 in. /ply) for a 
laminate thickness of 0.160 inch. The laminate configuration was [(±45)^, 90«, 

3 ^ ^ 

Ojq] . The density of the composite is about 0.054 pound per inch resulting in 
an airfoil weight of 18 pounds which is 2 pounds less than the design require- 
ment. The NASTRAN model is shown in figure 22 (a-suction surface and b-pressure 
surface) and consisted of 130 nodes and 120 CQUAD2 elements. The material prop- 
erties required for use in NASTRAN were generated using the composite mechanics 
code (ref. 1). Typical CPU time for the UNIVAC 1110 was: 7 minutes for the 

static case and 15 minutes for the first 10 vibration modes. 

NASTRAN predicted spanwise stresses on the suction surface are shown in 
stress contour plots in figure 23(a), and for the pressure surface in figure 
23(b). The maximum bending stress predicted by NASTRAN was about 7 ksi in the 
spanwise direction near the 1/3 chord point at midspan. The corresponding 
fracture stress measured in laboratory specimens with the same laminate config- 
uration was 128 ksi which is much more than five times the design load stress of 
7 ksi. Analogous comparisons were found for the other stresses. 

The first eight vibration modes predicted by NASTRAN are tabulated in ta- 
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ble V where the predominant vibration mode shape is also identified. The first 
predicted mode is 90 Hz which is 1.5 times the minimum design requirement of 
60 Hz. An interesting result from the vibration analysis is the second mode 
which is an accordion type mode shape. 

It is important to note that the laminate configuration described above was 
arrived at after two NASTRAN analysis cycles. This laminate configuration is 
not claimed to be optimum. Nevertheless, the example illustrates convincingly 
the effectiveness of NASTRAN as an analysis tool in the preliminary designs of 
composite airfoils. 


SUMMARY OF RESULTS AND CONCLUSIONS 

The major results and conclusions of this investigation are as follows: 

1. Out-of-plane eccentricities such as bending and twisting give rise to 
very large strains near the tab ends of off-axis composite tensile specimens 
and thereby initiate fracture at these locations. 

2. Vibration frequency measurements for evaluating the presence and/or ef- 
fects of small defects (about 2 percent by volume) in composites may not prove 
practical. 

3. Grip and/or end attachments can be designed to minimize transitional 
region strains due to load transfer in thin composite tubes using NASTRAN with 
plate elements. 

4. Effective bending moduli of composites with variable thickness can be 
determined using NASTRAN in conjunction with suitable frequency measurement ex- 
periments . 

5. NASTRAN can be used to help establish an effective plan for instrument- 
ing high-velocity- impact composite specimens. 

6. Stresses predicted using NASTRAN must be compared with in-situ material 
fracture stresses for reliable assessment of the design. 

7. The number of iteration analysis cycles in sizing airfoil-type composite 
components may be kept to as small as two with judicious use of NASTRAN. 

8. The results of the various analyses described in this summary paper dem- 
onstrate rather convincingly the versatility and effectiveness of NASTRAN as an 
analysis research tool in the rapidly advancing methodology associated with 
composite mechanics and composite structural components. 
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TABLE I. - FREE VIBRATION FREQUENCIES 

OF GRAPHITE /EPOXY (THORNEL 75 /epoxy 

[±45, 90^, 0„] ) COMPOSITE LAMINATE 
^ ^ s 


WITH DEFECTS 


Defect 

ratio 

(e) 

Frequency (Hz) for Mode 

1 

2 

3 

4 

0 

1385.9 

3907.0 

5453.9 

7819.4 

1/4 

1368.9 

3906.6 

5399.3 

7763.4 

1/2 

1321.5 

3903.5 

5345.0 

7589.1 

3/4 

1215.8 

3895.1 

5282.4 

7261.9 


a = defect width 
w = specimen width 


TABLE II. - CHANGE IN FREE VIBRATION FREQUENCIES 


CAUSED BY DEFECTS (THORNEL 75/EPOXY [±45, 9 O 2 , O^l^) 


Length 
of defect 

a 

Defect 

area 

(percent of 
total area) 

Percent change in frequency 
for mode shown 

1 

2 

3 

4 

1/4 

0.781 

- 1.2 

- 0.01 

-1.0 

-0.7 

1/2 

1.56 

- 4.6 

-0.09 

2.0 

-2.9 

3/4 

2.34 

- 12.3 

-0.3 

3.1 

-7.1 


a = defect width 
w = spedimen width 



TABLE III. - SUMMARY OF BENDING MODULI OF QUASI- ISOTROPIC 


COMPOSITES FROM GRAPHITE/EPOXY (MODMOR I/ERLA 4617) 


Layup 

Outer fiber 
direction 

Bending modulus* 
(106 psi) 

( 

o 

+1 

o 

0 

21.2 


90 

3.4 

(0, ±45, 90)^ 

0 

16.4 

s 

90 

3.4 

(0, ±30, ±60, 90) 

0 

17.6 

s 

90 

3.0 

(0, ±22.5, ±45, ±67.5, 90) 

0 

17.2 

s 

90 

2.8 


Average of two tests. 


TABLE IV. - COMPARISON OF COMPUTED 
AND FRACTURE STRESSES AT REGION 
OF MAXIMUM AVERAGE COMPOSITE 


STRESS AT DESIGN LOAD 


Blade point 
schematic 

Surface stresses, ksi 

Computed 

Fracture^ 

A(S.S.)^ 

(P.S.) 

-7.86 

67.1 



B(S.S.) 

(P.S.) 

-6.18 

64.16 


>.52.3 - 65.4 

C(S.S.) 

(P.S.) 

-4.75 

66.0 

J 



Measured in tensile specimens cut 
^ from fabricated blade. 

°For location of these points, see 
fig. 20; S.S. denotes suction 
surface; P.S. denotes pressure 
surface. 



TABLE V. - NASTRAN VIBRATION RESULTS FOR A FIBER 


COMPOSITE THIN SHELL AIRFOIL 


Mode 

ntamber 

Frequency , 
Hz 

Predominant mode shape 

1 

90 

First spanwlse bending (ISB) 

2 

118 

Accordian type 

3 

135 

Coupled torsion-bending (TB) 

4 

157 

Coupled spanwise-chordwise bending 

5 

181 

First torsion (IT) 

6 

187 

Second torsion (2T) 

7 

229 

Third torsion (3T) 

8 

245 

Second chordwise bending (2CB) 



0.057 in. 




^MiCARTA 

TAB 

Q/16in. 

THICK) 


r0.38in. 



Figure 1. ’ Specimen geometry and instrumentation of off-axis specimen i[0g]MOD I/E). 
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SPECliVlEN WIDTH. PERCENT 

Figure 4. - Out-of-plane bending effects on 
axial strain of the 30° off-axis specimen 
UOO in. -lb bending moment). 
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Figure 7. - NASTKAN 
model for composite 
specimens with de- 
fects; (297 nodes; 
256 elements). 


(a) NO DEFECT FREQUEN- 
CY: 7819 Hz. 


(b)WiTH DEFECT (a /w 
1/2) REQUENCY: 
7589 Hz. 


Figure 8. - Fourth free vibration mode shape of composite 
specimen with and without defect. 



Figure 9. - Photograph of instrumented tube mounted in grips, (MOD 1/E. [Og]). 
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Figure 11. - Schematic of 
NASTRAN model (Quadra* 
lateral plate elements). 
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Figure 14, - Typical composite cantilever specimen after high velocity Impact 
test (ref. 4). 
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Figure 14. - Typical composite cantilever specimen after high velocity impact 
test (ref. 4). 
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Figure 18. - Photograph of high-tip-speed composite blade. HTS/K601. (±40, ±20, 0). 


I 


I 



411 







X 



(3) SUCTION SURFACL 

Figure 22. - Airfoil NASTRAN model (130 nodes; 120 CQUAD2 elements). 
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(a) SUCTION SIDE. 

Figure 21 - NASTRAN predicted stress contour plots (ksi) in the 
spanwise direction on the outer surface of a composite airfoil 
(NACA 64A010). 




